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Matter Haloes in Filaments 



CO 

o 



o 

u 

Oh! 
O 



> 

O 

On 

in 

cn 
o 
m 



X 



M.A. Aragon-Calvo 1 * 

1 TTie Johns Hopkins University, 34-00 Charles St., Baltimore, MD, USA. 



March 2013 



ABSTRACT 

Dark matter haloes in cosmological filaments and walls have their spin vector 
aligned (in average) with their host structure. While haloes in walls are aligned with 
the plane of the wall independently of their mass, haloes in filaments present a mass 
dependent two-regime orientation. Here we show that the transition mass determining 
the change in the alignment regime (from parallel to perpendicular) depends on the 
hierarchical level in which the halo is located, reflecting the hierarchical nature of 
the Cosmic Web. By explicitly exposing this hierarchy we are able to identify the 
contributions of different components of the filament network to the spin alignment 
signal. We discuss a unifying picture to describe the alignment of haloes in filaments 
and walls consistent with previous results and our findings based on a two-phase 
angular momentum acquisition, first via tidal torquening and later via anisotropic 
mass accretion. 

The hierarchical identification and characterization of cosmic structures was done 
with a new implementation of the Multiscale Morphology Filter. The MMF-2 (in- 
troduced here) represents a significant improvement over its predecessor in terms of 
robustness and by explicitly describing the hierarchical relations between the elements 
of the Cosmic Web. 

Key words: Cosmology: large-scale structure of Universe; methods: data analysis, 
N-body simulations 



1 INTRODUCTION 

In the gravitational instability scenario galaxies acquire their 
angular momentum via tidal torque produced by a misalign- 
ment between the inertia tensor of the proto-halo and the 
embedding gravitational shear tensor generated by the sur- 
rounding matter distri bution. Thi s is the so c alled Tidal 
Torque Theory (TTT) (|Hovldll95ll : IWhitdll984l ). Although 
the TTT predicts a correlation between the spin orienta- 
tion of galaxies and their local Large Scale Structure ob- 
servations of spin alignme nts have been so far inconclu - 
sive and often co ntrad ic tory | Kashikawa & Okamura l (1 1992 
iNavarro etHI (120041); IPatiri et al.l (120061): iTrujillo et all 



On the theoretical/numerical front the picture is more 
clear. To date there is a consensus that the spin vector 
of haloes in walls lies in the plane of their host walls 
while haloes in filaments have their spin either paral- 
lel or perpen dicular to their ho s t filament depend i ng on 



their mass jAragon-Calvo et aL _ 2007 : 
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Libeskind et all 120131 ). An important aspect of the spin 
alignment of haloes in filaments is the existence of a tran- 
sition mass M tr , at which haloes change their orientation 
fro m parallel to perpen d icular . This effect was first reported 
by I Aragon-Calvo "e t al. (2007 ) and later confirmed by other 
authors dHahn et al.1 [20071 : IZhang et al.1 120091: ICodis et al.1 



, _ _ , .Hahn et aL 

I: ICodis et al.ll2012l: Trowlan d et al.1 



2007; 



2013; 



(2006 ): iLee Erdogdul (120071): iBrunino et al.1 
I Jones et al.1 (|2010l ): IVarela et al.1 (|2012l ): iTempel et al.1 



(|2013l ) among others). The reason is still not clear and may 
include the intrinsic difficulties in deducing the spin vector of 
galaxies, small galaxy samples, challenges in characterizing 
LSS from galaxy redshift catalogues, etc. 
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l2012l : ITrowland et"al](2013l : iLibeskind et alJl2013l ). The spin 
alignment of low-mass haloes also has a redshift depen- 
dence starting perpendicular and becoming pa rallel at z ~ 1 
(|Aragon-Calvo et al.ll2007l ; ITrowland et al.ll2013l ). This com- 
mon primordial orientation and subsequent change indicates 
an additional mechanism for angular momentum acquisition 
that has a stronger effect in low-mass haloes. In the obser- 
vational front the two regimes of alignment/anti-alignment 
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of low-mass and hi gh-mass haloes in filaments have been re- 
cently observed by Tem p el et al . (2013) if one assumes a cor- 
relation between Hubble type and galaxy/halo mass. They 
found that bright spirals have a weak tendency to have their 
spin parallel to their host filaments while Elliptical/SO show 
a strong perpendicular alignment. 

The existence of a hierarchy of structures in the Cos- 
mic Web in which larg e structures contain and are de - 
fined by sma l ler o nes (Sheth fc van de Wevgaertl (|2004h : 
Einasto et al.1 (|201lh ; Aragon-Calvo & Szalay (2013|) among 
others) begs the question of whether the halo alignment re- 
flects this hierarchy. The limited number of studies that link 
the alignment signal and the hierarchy of cosmic structures 
indicate that the transition mass of alignment in filam ents 
depends on the scale of the filament (jCodis et al. 2012) and 
the d istance from the halo to the closest cluster ([Zhang et al.l 
2009). However, to date no direct link has been stablished 
between halo alignment and the hierarchy of cosmic struc- 
tures. 

The TTT theory does not provide a mechanism (at least 
a simple one) to produce the observed two-regime alignment 
of haloes with their mass or to change the alignment of 
haloes from perpendicular to parallel with time. Recently an 
alternative mechanism for angular momentu m acquisition 
has b een discussed by Libeskin d et al ] (|2012h : ICodis' et al. 
(2012) and Trowland et al. I (|2013h in which angular momen- 
tum is generated by the anisotropic infall of matter. This 
is in term defined by the surrounding LSS morphology and 
dynamics, offering a natural process for spin-LSS correla- 
tions. This secondary-torque mechanism may help explain 
the change in alignment and possibly even contribute to the 
primordial torque generation. 



tures including dense haloes (below the smoothing scale). 
This reduces the dynamic range in the density field and 
greatly limits the contamination produced by dense haloes 
in the identification of filaments and walls. 

The identification of structures in the original MMF 
method and the one presented here is done based on the 
second-order local variations of the density field encoded 
in the Hessian matrix. LSS morphologies are associated 
to ratios between the eigenvalues of the H essian matrix 
(Ai < A 2 < A 3 , see lAragon-Calvo et all (|2007l )). The output 
of the MMF-2 is a set of "masks" sampled on a regular grid 
indicating which voxels belong to a given morphology at a 
given hierarchical level. For our purposes we generated a fil- 
ament, sub-filament wall and sub-wall masks (Fi, F2, Wi 
and W2 respectively, the subscript indicates the hierarchi- 
cal level). These masks roughly correspond to the "feature 
maps" T in the original MMF algorithm. A complete de- 
scription of the MMF-2 method will be presented in Aragon- 
Calvo (2013) in preparation. 



3 N-BODY SIMULATION AND DM HALOES 

The N-body simulations and related halo catalogues used in 
this w ork are based on the MIP correlated ensemble simula- 
tion. (| Ar agon- Cal vol [2012h . The MIP simulation consists of 
220 realizations sharing identical LSS fluctuations defined 
at k < k cut where k cut corresponds to a scale of 4 Mpc 
h _1 . At smaller scales (k > k cut ) realizations can be con- 
sidered independent between each other. Each realization 
in the ensemble contains 256 3 dark matter particles inside 
a box of 32 Mpc h _1 of side with cosmological parameters 
Q m 0.3, Q A 0.7, h 0.73, a 8 0.8. 



2 MMF-2: HIERARCHICAL IDENTIFICATION 
OF FILAMENTS AND WALLS 

In order to study the effect of environment in the alignment 
of haloes in filaments and walls we need to characterize the 
density field according to its local geometry, compute its lo- 
cal direction and, for the problem at hand, its hierarchical 
relations. The identification of filaments and walls was per- 
formed using a significantly improved implem entation of the 
Multiscale Morphology Filter (MMF, lAragon-Calvo et al.1 
<|2007h ). The MMF-2 performs the identification of struc- 
tures on a generalization of scale-space: the hierarchical- 
space, which exposes the hierarchical character of the den- 
sity field in contrast to the scale-space approa ch (as in the 
original MMF and its der ived implementations (|Zhang et al.l 
2009; ICautun et al1l2013h ) which emphasize the scale of the 
structures and is insensitive to their nesting relations. The 
hierarchical space used in the MMF describes both scale and 
nesting relations. 

We generate t h e hiera rchic al space as d e scribe d in 
lAragon-Calvo et al.1 (|2010bh and lAragon-Caivol (|2012h , by 
Gaussian-smoothing the initial conditions instead of the fi- 
nal density field as in the scale-space approach. This linear- 
regime smoothing is applied when all Fourier modes are in- 
dependent, allowing us to target specific scales in the density 
field before Fourier mode-mixing occurs. The smooth initial 
conditions evolve by gravity producing all the anisotropic 
features of the Cosmic Web but lacking of small-scale struc- 



3.1 Hierarchical space 

For the purpose of this paper we generated a two-level hi- 
erarchical space defined by linear-regime smoothing scales 
of 4 and 2 Mpc h _1 delineating the prominent struc- 
tures in the Cosmic Web and their substructure respec- 
tively. We denote these hierarchical levels as LEV-1 (fil- 
aments/walls) and LEV-2 (sub-filaments/sub- walls). From 
the templat e initial conditions us ed to generate the MIP en- 
semble (see lAragon- Calvol (|2012h for details) we generated 
low-resolution 128 3 initial conditions grids and smoothed 
them at the scales defined above. We then evolved the par- 
ticles to the present time z = using the PM-Tree code 
GADGET-2 (|Springell2005l ). These simulation were used only 
to compute density fields. 

3.2 Density fields 

From the particle distribution we c omputed densitie s fol- 
lowing the technique presented in lAbel et all (|2Q12h and 
IShandarin et al.l {2012). In this novel approach particles de- 
fine Lagrangian volumes that reflect the changes in den- 
sity as particles evolve deforming the original volume. La- 
grangian volumes are defined by 8 adjacent particles forming 
a "cube" that is divided into 6 tetrahedra. We estimate den- 
sities in the tessellation as the mean density of the adjacent 
tetrahedra to each vertex and interpolated to a regular grid 
taking care to account for multistreams. This method makes 
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Figure 1. The MMF-2 method. Top panels: hierarchical space. We show the volume rendering of the density field at z = for linear- 
regime smoothing at 4, 2 Mpc h _1 and no smoothing (left, center and rights panels respectively). The thickness of the slices was chosen 
to show as many structures as possible while avoiding confusion. Bottom panels: Hierarchical identification of filaments. The left and 
center panels show the density field inside filaments (Fi) and sub-filaments (F2) respectively. We highlight the filament mask as a 
semitransparent surface. The hierarchy of filaments (blue) and sub-filaments (red) is shown in the bottom right panel. 



density estimation computationally straightforward, it pro- 
duces a continuous volume-filling density field with practi- 
cally no artifacts in low and medium-density regions (our 
regions of interest) and is self-adaptive by construction. 

3.3 Halos 

Haloes were identified in each realization of the MIP en- 
semble using the friends-of-friends algorithm with a linking 
parameter b = 0.2. We imposed a count limit of 50 parti- 
cles corresponding to a halo mass of - 8 x 10 9 M h" 1 . We 
use the MIP ensemble in "stack mode", where halo cata- 
logues from all realizations are aggregated into one "master 
halo catalogue". The full MIP stacked ensemble contains 
— 5 x 10 6 haloes inside a volume of (32 Mpc h -1 ) 3 giving 
unprecedented halo number density. For each FoF group, 
we compute its mass, center of mass, mean velocity etc. The 
angular momentum was computed as: 



J = J^^i r i x (vi - v), 



(1) 



where the sum is over all the particles in the halo, mi is the 
particle mass, ri is the distance to each particle from the 
center of the halo, Vi is the peculiar velocity of the particle 



and v the mean velocity of the halo with respect to its center 
of mass. 



3.4 Spin alignment in filaments and walls 

The angle between the spin of the halo and its host filament 
or wall was computed as Of = 4>w and 0w — 90° — <pw, 
where <j)F —\ cos -1 (J • Ai) | and <fiw —\ cos -1 (J • A3) | are 
the angles between the spin of the halo and the direction 
of its host filament (Ai corresponding to the direction of 
constant density along the filament's axis) or the normal 
to the wall (A3 corresponding to the direction of maximum 
change in the density field). 



3.5 The filament /sub- filament sample 

Figure [T] shows the density field at z — from the two 
hierarchical levels defined by 4 and 2 Mpc h -1 linear-regime 
smoothing (LEV-1 and LEV-2 respectively) and the original 
initial conditions with no smoothing (top panels from left to 
right) . The original density field enclosed in the filament and 
sub-filament mask is shown at the bottom-left and center 
panels. The two-level hierarchy of filaments (the union of 
the two sets) is shown in the bottom-right panel. 
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Figure 2. Spin alignment of haloes as function of mass in filaments (left panel) and walls (right panel). The alignment is measured with 
the mean cos 9 where 0.5 corresponds to no preferred alignment. The sample was divided in haloes in filaments/walls (solid lines) and 
sub-filaments/sub- walls (dashed lines). Error bars per bin were computed from 1000 random realizations with the same number of points 
as their corresponding bin. 



In order to clearly differentiate between the alignment 
signal of haloes in filaments and sub-filaments we "cleaned" 
the sub-filament mask by removing the sub-filaments em- 
bedded inside filaments. We also constrained sub-filaments 
to be contained by walls. This removes the tenuous intra- 
void filaments. The mask of "clean" sub-filaments is then 
defined as: 

F 2 ,clean = F 2 D Wl D F? (2) 

where F 2 is the sub-filament mask, Wi is the wall mask 
and Ff is the complement of the filament mask. In what 
follows Fi and F 2 , dean denote the masks we used to com- 
pute spin-alignment of haloes in filaments and sub-filaments 
respectively 



4 RESULTS 

Our main findings are shown in Figure [2] where we mea- 
sure the mean cosine of the angle between the spin vector of 
haloes and the direction of their host filament/sub-filament 
(left panel) or wall/sub- wall (right panel) as a function of 
halo mass. A more detailed view of the spin alignment dis- 
tribution for haloes in filaments/sub-filaments is shown in 
Figure [3] Our results can be summarized as follows: 

• Haloes in filaments and sub-filaments present a two- 
regime spin alignment as a function of their mass: low-mass 
haloes are (in average) aligned parallel to their parent fila- 
ment while high-mass haloes are oriented perpendicular to 
it, confirming previous findings. 

• The mean alignment of low-mass haloes seems to con- 
verge at < cos# >^ 0.53 down to a halo mass of ~ 10 10 Mq 
h -1 in both filaments and sub-filaments. 

• The alignment of high-mass haloes does not seem to 



converge like in the case of the low-mass haloes. Instead, 
the alignment signal increases with increasing mass. 

• We found a significant difference in the transition mass 
M tr dividing parallel and perpendicular alignment for haloes 
in filaments and sub-filaments being M tr ,iow ^ 1.5 x 10 12 Mq 
h _1 and M tr , h i g h ~4x 10 11 M h -1 respectively. 

• The slope in the alignment-mass curve of high-mass 
haloes in sub-filaments is more pronounced than in the case 
of haloes in filaments. 

• Haloes in walls/sub- walls have their spin oriented in 
the plane of the wall for the full range of masses in our halo 
sample (~ 10 10 - 10 12 ' 5 M h" 1 ). 

• While the spin alignment distribution of haloes in fila- 
ments and sub-filaments seem scaled versions of each other 
(see Figure [3} the spin alignment distribution of low-mass 
haloes is practically the same for filaments and sub-filaments 
even though their transition masses are different by a factor 
of 3. 



5 CONCLUSIONS AND DISCUSSION 

Using an advanced LSS classification algorithm (MMF-2) 
and a novel N-body correlated ensemble simulation we were 
able to measure the halo spin alignment signal with unprece- 
dented detail in a variety of cosmic environments. Our re- 
sults confirm previous findings while significantly improving 
the alignment signal detection. However, we found a differ- 
ence in the transition mass M tr at which the spin alignment 
changes from being parallel (low-mass haloes) to perpendic- 
ular (high-mass haloes) depending on the hierarchical level 
of the host filament. This hierarchy-dependent effect has not 
been observed before due to the inability to explicitly target 
the nesting relations of the structures in the Cosmic Web. 
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Figure 3. Distribution of angles between the spin of a halo and 
its parent filament (black line) or sub- filament (red line). The 
samples were divided as high-mass haloes (solid line) and low- 
mass haloes (dashed line). The transition mass for filaments and 
sub-filaments corresponds to 6 x 10 11 and 2 x 10 12 Mq h _1 re- 
spectively. 



We found that haloes in walls have their spine aligned 
with the plane of the wall independently of the halo's mass. 

The explicit link between halo spin alignment and Cos- 
mic Web hierarchy was possible with the use of the MMF-2 
method (introduced here) that identifies cosmic structures 
and their hierarchical relations and benefits from an im- 
proved density estimation algorithm. The use of the MIP 
simulation allowed us to perform a high-resolution analysis 
of the density field over a relatively small volume while at 
the same time providing us with the large number of haloes 
needed to detect the spin alignment signal. 



5.1 On the origin and fate of the spin alignment 

The change of alignment of haloes in fila ments, the depen - 
dence of the transition mass with redshift ([Codis et al.| [2012) 
and filament hierarchy (this work) indicates a "secondary 
torque" process that acts cumulatively and affects all haloes 
but manifest itself first in low-mass haloes. Since tidal torque 
is mostly efficient before turnaround the process changing 
the spin orientation of haloes must be of a different nature. 
The anisotropic infall scenario offers a tantalizing alterna- 
tive. Walls and sub- walls accrete mass from their adjacent 
voids in the direction parallel to the normal of the wall in- 
ducing angular momentum on the plane of the wall so haloes 
retain their primordial spin alignment. On the other hand 
filaments and sub-filaments accrete mass mainly from their 
adjacent walls in the direction perpendicular to the filament 
thus inducing a parallel spin alignment. Low mass haloes are 
more affected because their low-mass makes then prone to 
have their spin direction changed for a given mass accretion 
rate compared to high-mass haloes. The stronger perpendic- 
ular alignment of haloes in sub-filaments in comparison to 
haloes in filaments reflects the more quiet environment char- 
acteristic of walls and their embedded sub-filaments. This is 
a natural consequence of the hierarchy of the Cosmic Web. 
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